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We review on the calculation of the heavy-quark photo-production vector form factors, with the full dependence 
on the mass of the heavy-quark. The Feynman diagrams are evaluated within the dimensional regularization 
scheme and expressed in Laurent series of e = (4 — D)/2, where D is the space-time dimension. The coefficients 
of the expansion are expressed in terms of harmonic polylogarithms. The numerical evaluation of the analytical 
formulas and the threshold limit of our results are presented. 



The measurement of the forward-backward 
asymmetry of the production of heavy quarks in 
e + e~ — > QQ processes, such as Ap B on the Z° 
peak 2], and, at the next generation of Linear 
Colliders, also the measurement of A l FB 0, is a 
stringent test of the Standard Model, providing a 
precise determination of sin 2 offl [3] . 

The present theoretical description includes the 
NNLO QCD corrections. For massless quarks, 
they were calculated numerically in 0] and an- 
alytically in 5 , while in [5] the order 0{ot 2 s )- 
corrections were calculated numerically for the 
b-quark retaining terms that do not vanish in 
the small-mass limit (constants and log-enhanced 
terms), but neglecting both terms containing lin- 
ear mass corrections, like m 2 /s, and terms in 
which such a ratio is enhanced by a power of 
log(s/m 2 ). In order to take into account also this 
kind of terms, a full analytic calculation in which 
the mass of the heavy quark is kept systematically 
different from zero is required. Three classes of 
contributions are involved: the tree level matrix 
elements for the decay of a vector boson into four 
partons, at least two of which being the heavy 
quark-antiquark pair; the one-loop corrected ma- 
trix elements for the decay of a vector boson into 
a heavy quark-antiquark pair plus a gluon; the 
two-loop corrections to the decay of a vector bo- 
son into a heavy quark-antiquark pair. 
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In this paper we consider the virtual 
NNLO QCD corrections to the process 7* — > 
Q(pi)Q{p2), with P = pi + P2 the momentum 
of the virtual photon and p\ , P2 the momenta of 
the outgoing on-shell quark and antiquark respec- 
tively (p 2 = p\ = — m 2 with m the mass of the 
heavy quark). The vertex amplitude V^{p\,p2) 
can be expressed in terms of two form factors as: 

V J, (pi,p 2 ) = u(pi)r^(p 1 ,p 2 )v(p 2 ), (1) 

r"(p llPa ) = Fi(pV + -LF 2 (p 2 ) a ^P„ (2) 
Zm 

where u(pi),v(p2) are the spinor wave functions 
of the quark and antiquark, = — §[7^,7"], 
and p 2 = P 2 /m 2 = —s/m 2 , with s the squared 
cm. energy. 

At the two-loop level, the contributions to 
Fi(p 2 ) and Fi{p 2 } come from the diagrams shown 
in Fig. n] They are evaluated in jZj using a 
method based essentially on two steps: 

1. By means of suitable projector operators, 
the contributions to F\{p 2 ) and F2{p 2 ) 
are expressed in terms of a big number 
of scalar integrals, whose ultraviolet (UV) 
and infrared (IR) singularities are regular- 
ized within the dimensional regularization 
scheme. Applying the so-called Laporta al- 
gorithm, introduced in [H] and become, in 
the last few years, a standard technique for 
the calculation of Feynman diagrams, all 
these scalar integrals are expressed in terms 
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Figure 1. Two-loop Feynman diagrams. The 
curly lines represent gluons; the double straight 
lines, quarks of mass m; the single straight lines, 
massless quarks and the dashed lines ghosts. 



of 34 independent integrals called Master 
Integrals (Mis). The algorithm is based on 
Integration by Parts Identities UJ, Lorenz 
invariance identities |10j and General Sym- 
metry Relations. It is completely algebraic 
and it is implemented in a code written in 
FORM [TT] 

2. The calculation of the Master Integrals, 
then, is carried out (12113) by means of the 
Differential Equations technique |14j . 

The renormalization of UV divergences is per- 
formed in a hybrid scheme: we renormalize the 
mass and wave function of the heavy quark in the 
on-shell renormalization scheme, while the cou- 
pling constant, as well as the light-fermion and 
gluon wave functions, are renormalized in the MS 
scheme. 

The final result of the calculation is an ana- 
lytical expression for the form factors, written in 
terms of harmonic polylogarithms of one variable 
|15I16| . The result has still IR divergences that 
appear as poles in 1/e. These divergences have 
to be canceled against those arising from the real 
radiation, which are not taken into account here. 

We present, in this short review, the numerical 
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Figure 2. One-loop form factors, F± and F2 , in the 
case of top-antitop production. The charge form 
factor (upper plot) has an IR pole, while the mag- 
netic form factor is IR finite. Note the Coulomb 
divergences of the form factors in the threshold 
limit . The real part of the pole of F\ and the real 
part of F 2 are, instead, finite (0 and -Cp respec- 
tively). 

evaluation of the different parts of the form fac- 
tors (real and imaginary parts of the poles and 
finite parts) at one- and two-loop level above the 
threshold for the case of the production of a pair 
of top-antitop quarks (one massive and five light 
flavours). The numerical evaluation is performed 
in F77 using the subroutine of ^Hl- The color 
charges are set according to the SU(3) group. 
The cm. energy runs from the threshold value 
(y/s ~ 2m, m = 178 GeV) to 2 TeV. The plots 
are shown in Figs. |21 El and El 

The knowledge of the exact analytical formula 
for the form factors allows a perfect control on 
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Figure 3. Two-loop charge form factor, F\, in the 
case of top-antitop production. The charge form 
factor has a double IR pole. 
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Figure 4. Two-loop magnetic form factor F2, in 
the case of top-antitop production. The magnetic 
form factor has a simple IR pole. 



their numerical evaluation as well as on their an- 
alytical behaviour in all the regions of the spec- 
trum. As an example, consider the threshold re- 
gion. The NNLO determination of the thresh- 
old limit of the vector form factors were already 
present in the literature J7| as an independent 
calculation, performed with a different technique. 
With the method outlined above, however, this 
result can be achieved as a limiting case of the 
full calculation. 

In fact, let us define the real and imaginary 
parts of the UV-renormalized form factors at one 
and two loops as: 



F, 



(1L21) 



C(e) ^f\ 11 



21) 



in^sF, 



(11,21) 



}, (3) 



with £ = 1,2 and C(e) = (47r) e r(l + e). More- 
over, consider the total cross-section for the pro- 
cess e + e~ — * QQ near the threshold: 



*QQ Z 



I- 



2tt \2ttJ 



(4) 



where cr (0) = {2ira 2 e 2 Q /s) /3(3 - (3 2 ) is the tree- 



level cross-section and /3 = y 1 — is the quark 

velocity. The one- and two-loop corrections, A^ 1 ) 
and A( 2 \ can be written in terms of the form 
factors, up to 0(/3 2 ), in the following way: 



A (2) 



(i0 



(5) 



^f\ 



2$tF[ ll) UF^ l) 
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2^F\ 



(21) 
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QF. 



(ii) 



(0) 



Using the analytical expressions for the threshold 



limits of ft^J"', QF^ 1 ', Ztf^' 1 , <Sf£ 1) , nFl 
and $tF^ l \ given in |7], the functions A^ and 
A< 2 ) read: 

A W=C F [^-8j, (7) 
AW = C 2 |±S-^ + 39 



7 (2l) 
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151 
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32 
T 



C 2 -8C 2 ln(2/3) 

C2 



44 



(8) 



where £ 2 and £3 are the Riemann C, function eval- 
uated in 2 and 3 respectively. Eqs. (|7I8|1 are in 
complete agreement with the results in |17| . 

To conclude, we presented the calculation of 
the two-loop QCD contributions for the vector 
form factors of the photo-production of heavy 
quarks. By means of the Laporta and differential 
equations methods we were able to calculate ana- 
lytically the form factors, keeping the full depen- 
dence on the heavy-quark mass. We presented the 
numerical evaluation of the form factors as well as 
their threshold limit, interesting for the calcula- 
tion, in this kinematical region, of the total cross 
section for the photo-production of heavy quarks. 
This result is in complete agreement with the ex- 
isting result in the literature. 
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